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ABSTRACT 


A  noise-generating  algorithm  and  associated  computer  program 
for  well-defined  testing  of  beamformers  are  described*  The  algorithm  is 
especially  suitable  for  superdirective  arrays  of  underwater  hydrophones 
as  it  generates  Gaussian  noise  of  specified  coherency.  Statistical 
properties  of  the  generator  are  confirmed  to  be  those  planned,  and  the 
ability  of  the  generator  to  synthesize  noise  for  isotropic  or  surface 
noise  sources  is  verified  for  three-element  arrays.  Cumulative 
distributions  for  estimated  coherency  were  obtained  for  the  model. 
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INTRODUCTION 


Computer  programs  for  theoretical  testing  and  comparison  of 
beamforraing  algorithms  require  noise  generating  algorithms  that 
synthesize  noise  of  known  coherency  and  statistical  properties. 

There  is  a  significant  advantage  in  using  noise  synthesizers 
to  select  suitable  beamformers  economically  before  field  testing.  The 
type  of  noise  generated  can  be  controlled  and  the  beamformers  tested  for 
a  set  of  defined  and  reproducible  noise  conditions.  A  considerable  time¬ 
saving  results  since  the  testing  of  the  beamformers  for  noise  conditions 
that  might  be  met  in  the  field  over  several  years  can  be  done  in  the 
laboratory  in  a  matter  of  days. 

For  arrays  of  widely  spaced  sensors,  where  the  noise  is 
uncorrelated  from  sensor  to  sensor,  noise  generators  simply  consist  of 
uncorrelated  noise  sources,  one  noise  source  for  each  sensor.  However, 
for  arrays  of  closely  spaced  sensors,  a  model  to  generate  noise 
correlated  from  sensor  to  sensor  is  required.  This  memorandum  describes 
the  simulator,  verifies  its  statistical  properties,  and  delineates  those 
noise  fields  that  can  be  represented  by  the  simulator. 

THEORY 


A  beamformer  that  explicitly  includes  a  device  to  calculate 
Fourier  transforms  of  the  hydrophone  outputs  is  shown  in  Figure  1.  For 
computational  efficiency,  the  noise  generator  described  here  produces  the 
Fourier  transforms  of  the  noise  directly,  instead  of  generating  the  time 
series  of  the  noise  and  subsequently  calculating  the  transform.  These 
transforms  are  arranged  to  be  random  variables  with  a  Gaussian 
distribution  that  has  been  found  to  be  characteristic  of  ambient  noise 
over  intervals  of  a  few  minutes** 


To  generate  noise  of  specified  coherencies  between  the  n 
sensors,  the  Fourier  transform  Xj^w),  of  the  ith  sensor  at  the  frequency 
a) ,  is  written  as  a  linear  combination  of  real  and  imaginary  pairs  of 
Gaussian  distributed  random  variables  Z^u).  Both  the  real  and  imaginary 
parts  have  a  mean  of  0  and  a  variance  of  0.5.  Dropping  reference  to 
frequency,  these  linear  combinations  are  written: 


X1  =  all  Z1  +  a12  z2  + 


+  aln  zn 


X2  -  a2i  +  a22  Z2  + 


+  a2n  zn 


(1) 


xi  -  ail  Z1  +  ai2  z2  + . +  aln  Zn 
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Figure  1.  In  the  generalized  beamforraer  shown  the  time  series  x^t)  is 
Fourier  transformed  to  X^(  to  j)  and  the  transforms  are 
multiplied  by  the  weights  j 


The  values  of  the  a^j,  which  are  restricted  to  be  real,  are 
determined  by  the  requirements  that  on  the  average  the  noise  field  power, 
qij  i*»j,  be  homogeneous  (the  same  at  all  hydrophones  and  equal  to  unity) 
and  that  the  average  noise  field  coherency,  q^j  i/tj,  between  sensor  pairs 
be  as  specified  by  the  user  (e.g.  isotropic  noise).  These  two  conditions 
may  be  written 

^ij  *  ^i^j  i , j*l , 2 , . • .n.  (2) 

In  addition,  the  simplifying  assumption  was  made  that 

atj  -  0  j>i.  (3) 

By  combining  (1),  (2),  and  (3)  and  using  the  independence  of 
the  it  can  be  shown  that 

_  i 

q  jj  *  XiXj  *  I  aika  jk  j*  1 , . .  i  f  i^"  1 ,  •  •  *n •  ) 

k*l 

These  equations  are  solved  for  a^j  and  the  Fourier  transforms 
are  then  calculated  from  Equation  (1).  A  listing  of  the  noise 
generating  program  is  contained  in  Appendix  A.  The  subroutine  Gauss  4 
called  by  the  noise  generator  has  been  extensively  tested  and  found  to  be 
faster  computationally  and  better  statistically  than  the  random  number 
generator  'Gauss*  supplied  with  IBM  systems  software^. 

The  noise  generating  algorithm  cannot  solve  for  a^j  for  all 
arbitrary  sets  of  coherency  values.  Firstly,  the  form  of  Equation  (3) 
restricts  noise  fields  modelled  to  those  for  which  Qij^ji*  By  doubling 
the  number  of  random  variables  ,  complex  qAj  could  be  accommodated. 
Secondly,  even  for  a  three-element  array  the  requirement  that  333  be  real 
restricts  permissible  q^j.  To  obtain  some  indication  of  whether  this  is 
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a  severe  limitation,  examples  of  noise  fields  that  give  real  333  for  a 
three-element  'equispaced1  horizontal  line  array  were  determined 
numerically  and  theoretically. 

The  condition  on  q^j  that  must  be  satisfied  for  real  333  for 
any  three-element  array  is, 

q13  q23  +  2q13  q23  q12  +  qi2  “  1  -  0  (5) 

This  condition  is  a  special  case  of  the  more  general  requirement  that  the 
cross  spectral  matrix  be  Hermitian  positive  semidef inite^*  Equation  (5), 
which  is  derived  in  Appendix  B,  was  tested  for  isotropic  noise,  i.e. 
noise  whose  coherency  is  given  by 

slnUdtj)  _  (6) 

kdjj 

and  for  surface-generated  noise  for  which  the  coherency  can  be  expressed 
as 


q  =  2lDm!  Vkd> 

J  (kdij)” 


(7) 


where  k  is  the  wave  number,  djj  is  the  sensor  separation,  and  Jm  is  the 

Bessel  function  of  the  first  kind  of  order  m.  The  condition  specified  by 

Equation  (5)  is  satisfied  for  three-element  equispaced  arrays  for 

isotropic  noise  and  for  surface  generated  noise  for  m  =  0,  1,  2  and  —  up 

X 

to  0.95.  This  was  shown  theoretically  for  surface  noise  as  outlined  in 
Appendix  C  and  numerically  for  isotropic  noise.  Beyond  0.95  of  a 
wavelength  the  model  approaches  that  of  independent  noise  sources,  one 
noise  source  for  each  hydrophone. 

It  might  be  thought  that  allowing  a^j  to  be  complex  would 
remove  the  restriction  imposed  by  Equation  (5)  and  allow  modelling  of  a 
wider  range  of  noise  fields.  However,  even  for  complex  a  ^  j  the 


restriction  on  the  noise  coherency  as  defined  by  Equation  (5)  remains* 
Furthermore,  allowing  a^j  to  be  complex  introduces  a  new  difficulty. 
While  for  real  a^ j  all  sensors  will  have  a  uniform  distribution  of  the 
phase  shift  between  the  real  and  imaginary  parts  of  the  Fourier 
transform,  complex  a^j  introduces  the  situation  where  there  are 
distinctly  different  distributions  for  different  hydrophones;  this  is 
equivalent  to  saying  that  the  noise  field  is  not  homogeneous  in  the  phase 
shift  distribution  and  is  therefore  rather  unrealistic.  The  restriction 
to  real  is  thus  not  purely  arbitrary. 

D  CUSSION  OF  RESULTS 

Tests  were  carried  out  to  determine  whether  the  synthesizer 
produced  noise  with  the  desired  statistical  properties.  Firstly,  the 
Kolmogorov-Smirnov  test  was  applied  to  test  the  hypothesis  that  the 
Fourier  transform  amplitudes  are  Gaussian  distributed  random  variables. 
The  test  was  applied  to  the  cumulative  distribution.  Each  cumulative 
distribution  tested  contained  500  samples  of  the  transform  and  100 
cumulative  distributions  were  tested.  A  significance  level  was 
calculated  for  each  of  the  100  cumulative  distributions.  The 
significance  level  indicates  the  probability  that  the  cumulative 
distribution  would  have  occurred  by  chance.  Individual  significance 
levels  were  consistent  with  the  hypothesis  that  the  sample  came  from  a 
population  of  Gaussian  distributions. 

The  100  significance  levels  from  the  Kolmogorov-Smirnov  test 
were  also  examined.  They  lie  between  0  and  100%  and  should  have  an  equal 
probability  of  occurrence,  i.e.  the  significance  levels  should  be 
uniformly  distributed.  The  observed  set  of  100  significance  levels 
obtained  in  the  Kolmogorov-Smirnov  test  departed  somewhat  from  a  uniform 
distribution.  It  was  necessary  to  know  whether  this  departure  from  a 
uniform  distribution  was  likely  to  occur  by  chance.  Again  the 
Kolmogorov-Smirnov  test  was  used  to  investigate  the  hypothesis  that  the 
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significance  levels  were  uniformly  distributed.  This  hypothesis  of 
uniform  distribution  could  not  be  rejected  at  the  27%  level,  i.e.  there 
is  approximately  one  chance  in  four  of  obtaining  this  particular 
distribution  or  one  with  a  greater  deviation  from  uniformity.  Thus  there 
is  no  reason  to  suspect  the  original  hypothesis  of  the  Fourier  transform 
amplitudes  being  Gaussian  distributed.  Indeed  confidence  in  the 
hypothesis  is  increased. 

Secondly,  the  power  from  each  sensor  was  tested  to  determine 
whether  the  power  was  chi-squared  distributed  with  two  degrees  of 
freedom.  Significance  levels  were  calculated  from  the  Kolmogorov-Smirnov 
test  for  cumulative  distributions  containing  100  samples  of  the  power  in 
20  trials  with  5  sensors.  The  calculated  individual  significance  levels 
were  consistent  with  the  chi-squared  hypothesis.  Again  to  aid  in  the 
evaluation  of  the  significance  levels  as  a  group,  the  hypothesis  that  the 
significance  levels  were  uniformly  distributed,  as  they  should  be,  was 
tested  with  the  Kolmogorov-Smirnov  test.  It  was  found  that  the 
hypothesis  could  not  be  rejected  at  the  77%  level.  These  results  are 
taken  as  confirmation  that  the  power  is  indeed  chi-squared  distributed 
with  two  degrees  of  freedom  as  was  intended. 

Thirdly,  the  phase  angle  of  the  sensor  outputs  should  be 
uniformly  distributed.  In  the  20  trials  with  5  sensors,  significance 
levels  were  calculated  using  the  Kolmogorov-Smirnov  test  for  cumulative 
distributions  containing  100  samples  of  the  phase  angle.  Again  the 
individual  significance  levels  were  consistent  with  the  hypothesis  under 
test.  Since  the  significance  levels  should  themselves  be  uniformly 
distributed,  they  were  tested  for  a  uniform  distribution  with  the 
Kolmogorov-Smirnov  test.  The  hypothesis  of  a  uniform  distribution  of  the 
significance  levels  could  not  be  rejected  at  the  97%  level  so  that  the 
hypothesis  that  the  phase  of  the  sensor  output  is  uniformly  distributed 
gains  further  support. 


PROBABILITY 


Additional  checks  were  made  to  verify  that  the  algorithm 
produced  noise  whose  coherencies  converged  to  the  specified  coherence  for 
the  noise  field*  Hydrophone  outputs  were  synthesized  for  isotropic  noise 
and  also  for  a  surface  noise  field  represented  by  JQ(kd)  as  given  by 
Equation  (7)  for  m^O.  This  was  carried  out  for  up  to  five  hydrophones 
for  various  sensor  configurations  and  in  all  cases  solutions  were  found 
for  the  aij*  The  calculated  coherencies  for  estimates  made  from  samples 
of  100  coherencies  produced  by  the  simulator  showed  a  bias.  That  bias 
agreed  well  with  the  bias  given  by  Benignus^  for  coherencies  generated 
from  two  independent  Gaussian  noise  sources. 

Cumulative  distributions  for  the  coherencies  were  calculated 
for  a  sample  size  of  100  at  9  selected  coherencies.  These  are  plotted  in 
Figure  2  to  characterize  the  model  and  enable  comparison  of  measured 
cumulative  distributions  of  coherency  with  coherency  calculated  from  the 
model.  For  sample  sizes  between  2  and  100  the  95%  confidence  limits  are 
summarized  in  Figure  3* 


SQUARED  COHERENCY 

Figure  2.  Cumulative  frequency  distributions  for  the  calculated  mean 
squared  coherency.  To  obtain  the  curves  plotted,  500 

estimates  of  coherency  were  made  with  a  sample  size  of  100* 
The  true  squared  coherency  is  listed  beside  each  curve. 
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CONCLUSIONS 


The  algorithm  meets  the  requirement  of  generating  noise  for 
testing  beamformers  for  closely  spaced  arrays.  This  enables  testing  and 
comparison  of  beamformers  in  the  laboratory  for  noise  fields  of  defined 
and  reproducible  properties. 

It  was  verified,  for  three-element  equispaced  arrays,  that 
the  algorithm  is  able  to  model  noise  fields  with  coherencies 
corresponding  to  isotropic  noise  and  to  surface  noise  fields.  However, 
the  algorithm  does  not  generate  noise  for  all  arbitrary  noise  fields#  An 
expression  that  must  be  satisfied  by  the  coherencies  for  a  three-element 
array  was  obtained. 

The  statistical  properties  of  the  synthesizer  were  confirmed 
to  be  those  for  Gaussian  noise  and  cumulative  distributions  of  the 
coherency  were  obtained. 
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APPENDIX  B 


In  this  appendix  the  condition  on  the  noise  coherencies  qAj 
for  real  .133  Is  derived  for  a  three-element  array.  As  previously  the 
hydrophone  output  XA  is  written 
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APPENDIX  C 


In  this  appendix  some  noise  fields  that  nan  be  modelled  by 
the  algorithm  are  determined*  The  investigation  is  limited  to 
three-element  'equispaced1  horizontal  arrays*  For  an  equispaced  array 
qi2  ■  <123  an<*  becomes, 

q13  "  1  ~  2q12  ^q  13  ^  -  0 

for  333  real.  This  equation  may  be  written 

(q13  -  l)(q^3  +1-2  q22)  £  0 

and  since  (^13  -  1)  is  always  negative  we  require 

2  q12  "  q13  '  1-°  (C1) 


for  real  333. 

Case  1 

For  surface  noise  whose  coherency  ran  be  represented  by  JQ(x) 
where  x  *  kd,  the  left-hand  side  of  (Cl)  becomes 

2J02(x)  -  J0(2x)  -  1  (C2) 

To  evaluate  this  expression  we  have  the  addition  theorems  for  Bessel 
function*^: 

<x> 

J  2(X)  +  2  I  J2(x)  -  1  (C3) 

k-1  R 

00 

and  J0(2x)  »  JQ2(x)  +  2  E  (-I)*  Jjl2(x)  (C4) 

k»l 


02 


Substituting  (C4)  In  (C2)  and  splitting  the  sum  into  even  and  odd  parts 
we  obtain 
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and  by  applying  (C3) 
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J2k(x) 


This  verifies  that  the  left-hand  side  of  (C2)  is  certainly  less  than  or 
equal  to  zero  for  all  x.  Thus  the  algorithm  can  find  real  aj3  and 
synthesize  acoustic  noise  for  surface  noise  of  the  form  JQ(x)  for  all 
hydrophone  separations  with  a  three-element  equispaced  array. 


Case  11 


For  surface  generated  noise  fields  the  noise  coherency  can  be 
expressed  by^ : 


qij  - 


2“m!  Jm(x) 
x® - 


oo 
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To  simplify  substitution  into  (Cl),  the  test  for  real  333,  we  note  that 
|J*(x) 


q12 


x2  ml  V*)  +  +  (~Ok  x2km!  Jm(x) 
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-  4x2ml  + . +  4(-l)k  x2k  ml  +  ... 

22<“  +1>!  22k  k!(m  +  k)! 


(C7) 


03 


Now  substituting  in  (Cl),  grouping  even  and  odd  terms  and  using  £  to 
denote  the  even  terms,  the  left-hand  side  of  (Cl)  becomes 


are  negative  for  all  x.  The  third  term  is  negative  provided  x  <  6*  This 
implies  that  333  is  known  to  be  real  under  the  following  conditions. 


1.  the  array  consists  of  three  equispaced  sensors  in  a  line; 

2.  the  noise  field  is  of  the  form  (C5); 

3.  the  largest  hydrophone  separations  are  £  0.95 

wavelengths. 

It  was  also  found  from  numerical  evaluation  of  Equation  (Cl)  that  333  is 
real  out  to  hydrophone  separations  of  1.5  wavelengths  for  m  -  1,  2,  or  3 
with  surface  noise  fields  of  the  form  given  by  (C5). 
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